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Abstract 

The availability of controllable macroscopic devices, which maintain quantum coherence over 
relatively long time intervals, for the first time allows an experimental realization of many effects 
previously considered only as Gedankenexperiments, such as the operation of quantum heat engines. 
The theoretical efficiency n of quantum heat engines is restricted by the same Carnot boundary rjc 
as for the classical ones: any deviations from quasistatic evolution suppressing n below r\c- Here 
we investigate an implementation of an analog of the Otto cycle in a tunable quantum coherent 
circuit and show that the specific source of inefficiency is the quantum squeezing of the thermal 
state due to the finite speed of compression/expansion of the system. 
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I. QUASISTATIC QUANTUM OTTO CYCLE 



Quantum heat engines (QHE) are engines with a quantum coherent working body. After 



they were introduced Q, l2| as the generalization of classical heat engines for lasing, they 
were extensively used in Gedankenexperiments helping clarify subtle points of statistical 



mechanics, like the ro 
regime (see, e.g 



1C 1U 



e of fluctuations and the operation of Maxwell's demon in the quantum 
). Even though QHEs demonstrate some very nonclassical features, of 
course these do not violate the second law of thermodynamics. In particular, QHEs satisfy 
the Carnot inequality, 

R T 

(1) 



V = -Vc 



T 

1-— , 

Th 



where R is the mechanical work performed per cycle, Qh is the energy obtained from the 
heater (an equilibrium reservoir at temperature Th)} and Q c is the energy transferred to the 
cooler (at temperature T c ). 
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13), 



The latest developments in solid state-based quantum computing (see, e.g., 
especially with superconducting devices, gave this line of investigation more experimental 
relevance. Unlike the various mesoscopic cooling schemes in operation since the mid-1990s 



see [la]), the sideband cooling of superconducting qubits, achieved in [19 



does not 



involve the exchange of particles between the system and the thermal reservoirs and is 
therefore much closer to Gedankenexperiments' QHEs. 

Any heat engine operating along any other thermodynamic cycle is less efficient than the 
ideal Carnot engine. Moreover, any heat engine operating at finite speed is less efficient than 
the same engine operating quasistatically. The source of inefficiency is in the irreversibility of 
non-quasistatic processes, but the specific mechanism for the QHE depends on the particular 
realization and is often assumed to be due to the loss of quantum coherence. In this paper 
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l25|) and show that the source 



we consider the quantum Otto cycle (see, e.g., 
of its inefficiency at finite operation speed (compared to the quasistatic case) lies in the 
reversible, quantum coherent squeezing of the quantum state of the working body. 

Let us consider a textbook situation, where there are two thermal reservoirs at tempera- 
tures T c , Th, and a system which can be periodically put in thermal equilibrium with either 
of them (e.g., 26j). In its quantum analog, the reservoirs and the system are represented 
by three harmonic oscillators, with frequencies 0J c ,Uh, and u, respectively, assuming that uj 
can be changed at will. In the following, the nature of these oscillators is immaterial. For 
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the sake of definitiveness, we will talk about superconducting LC-circuits with a Josephson 
junction, which can be biased by a current or an external magnetic flux. The tunability of 
such a system in the quantum coherent regime was repeatedly demonstrated in experiments 



[e.g., 



n 



28|). We assume that 

u h > u c , (2) 



29| . the occupation 



but will not impose any restrictions on the baths temperatures. In such an implementation, 
the mechanical work produced by, or applied to, the system is electrical and is determined 
by the quantum state-dependent energy required to change the bias current in a Josephson 
junction. 

The quantum Otto cycle runs as follows (Fig. [I]). The system is initially tuned to and 
equilibrated with the "cold" reservoir (A). Then it is quasistatically tuned in resonance 
with the "hot" reservoir. According to the quantum adiabatic theorem 
numbers of its energy eigenstates will not change during this process, while mechanical 
work must be performed to increase u. At point (B) the system is in resonance, but out 
of equilibrium, with the "hot" reservoir, u = Uh- After the equilibration is achieved (C), 
the system is quasistatically brought back in resonance with the "cold" reservoir, and is 
again equilibrated with it at (A), closing the cycle. On the adiabatic stages AB and CD, 
the system is isolated from the outside world (neglecting the finite, but ideally vanishing 
linewidth of the system and reservoirs' energy levels), and all the energy exchange takes 
place along the "isochoric" stages BC and DA. The work performed by the system is given 
by the area enclosed by the contour ABCD: 

rB rD 

R= (n) du - / (n) duj = ±Aabcd, (3) 

J A JC 

where (n) is the expectation value of the photon number in the system. The role of volume 
is played by the inverse frequency, 1/u, and the pressure in the system is (n)u> 2 . For the 
clockwise sense R > 0, and the device works as a heat engine; otherwise R < 0, and it is a 
refrigerator, transferring energy from the "cold" to the "hot" reservoir. 

The calculations for the ideal Otto cycle are straightforward: denoting by \n,c(h)) the 
nth energy eigenstate of an oscillator with frequency u c (uh), and by Z c ^ the corresponding 
partition function, we can write for the density matrix of the system at points A, B, C, D 

Pa = Z~ l ^2e~^\n 1 c)(n ) c\ 1 p B = Z~ x ^e""^|n, h)(n, h\, 
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FIG. 1: (Color online.) The expectation value of the photon number (n) in the tunable oscillator 
("working body") versus the oscillator frequency uj. This quantum Otto cycle is for a system com- 
prising of two oscillators in thermal equilibrium ("reservoirs") and a tunable oscillator ("working 
body"). This cycle is realized by changing the latter 's frequency between oj c and ujh along AB and 
CD (adiabatic stages). Along BC and DA (isochoric stages), the tunable oscillator equilibrates with 
the reservoirs. Whether the system operates as a heat engine (red, upper cycle) or a heat pump 
(blue, lower cycle) is determined by the parameter A = u c /T c — lu^/T^, Eq. (J7J). (a) Quasistatic 
case, (b) Finite speed case [calculated using Eq. (|2Q|) with the oscillator frequency w as a function 
of time shown in the inset of Fig. [5]. 

Pc = Z^ 1 J2 e ~~\n,h)(n,h\, p D = Z^ 1 ^ e~~ \n, c) (n, c\ , (4) 

n n 

which immediately yields the expressions for the energy of the system at the corresponding 
points and for the work done and the energy received by the system: 
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n 



(6) 



R\ = E A — E B = —(u h — lo c )u c , R 2 = E c — E D = (u h — io c )v h , 

R = Ri + R 2 = (u h - w c )(v h - p c ); (5) 
Qt = E c - E B = u h (i> h - u c ); Q 2 = E A - E D = -u c (u h - u c ), 

where Vi = n(x«) + | (i = c,h), and the thermal population of the oscillator's energy levels, 
fi(xi), has the standard form: 

'uA = 1 
.Tj/ exp [Ui/Ti] - 1' 

Below we will also use k = (2z/) _1 = tanh(w/2T). 

As expected, we see from Eq.([6]), that here always -R1-R2 < and Q1Q2 < 0. The direction 
of the cycle is determined by the parameter 

A = 7f - t£- (7) 

-Lc J-h 

For A > 0, the net work R = Ri + R 2 > 0, and the device works as a heat engine with 
efficiency 

v = — = l-—. (8) 

VI 

For A < 0, it becomes a refrigerator with efficiency 

C = ^ = ^— (9) 

It is also straightforward to check that the Clausius inequality, C = AS — AQ/T > 0, is 
satisfied on the non-adiabatic stages of the cycle. Indeed, the entropy on the compression 
stage is 

Sab = In (v c - -) + ^ (v e + -) = - ln[2 sinh(w c /2T c )] + (u c /2T c ) coth(u; c /2T c ), (10) 

and for that on the expansion stage, Sqd, we replace (u c /T c , u c ) with (uh/T h , u h ). From 
here, using the identity 



Uk-u e = 7j v (11) 

sinh (*h) sinh (^-) 



and, for positive x, 

ln(sinha;/ sinh?/) < (x — y) cothy, 



we find that indeed 

Cbc — Sqd — Sab — t)t~ — 0; 

-Lh 

Cda = Sab — Sqd — t)t~ > 0. (12) 

J- c 

An interesting feature of the quantum Otto cycle is that, unlike the Carnot cycle, its 
efficiency is independent on the temperatures of the heater and cooler (cf. Q). The Carnot 
inequality is not violated though, since precisely at the point where r\ = rjc, the parameter 
A goes through zero and switches sign, turning the heat engine into a refrigerator. However, 
the independence of the cooling efficiency ( on temperature may make a quantum Otto 
fridge a candidate for the realization of self-cooling qubits. The minimal temperature such 
a fridge can cool down the "cooler" is limited by the condition A < 1 to 

T c = T h — . (13) 

I — I 

The minimal temperature that can be achieved using other methods is discussed in [301 ] . 



II. QUANTUM OTTO CYCLE AT FINITE OPERATION SPEED 



In order to analyze the effects of a finite speed of operation of the Otto engine, it is con- 
venient to use the Wigner function in the basis of coherent states of the harmonic oscillator 
(e.g., 



3l|), which was used to investigate parametric squeezing 
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32J. It is well known 



that a sudden change of frequency of a linear oscillator transforms a coherent state into a 
squeezed state, with the squeezing proportional to the ratio of initial and final frequencies, 
while a quasistatic change does not 0-0. 



Note that squeezing is a unitary ( i.e., reversible ) operation. The von Neumann entropy 
of a quantum state, S[p] = — trplnp, is not changed by squeezing. The situation is different 
for the so-called "energy entropy" (see, e.g., [22j), 



S f 



J2Pn lll Pn > S[p], 



(14) 



n=0 



where p n are the diagonal components of the density matrix in the energy representation. 
Note that Se and S coincide only if [p, H] = 0. The decoherence processes, which eventually 
eliminate the off-diagonal elements of the density matrix in the energy basis, thus increase the 
von Neumann entropy and can be associated with the inner friction processes. The energy 



entropy is therefore a quantitative measure of the inner friction in a quantum system, which 



is appropriate for the analysis of quantum heat engines operating at finite speed [22] . The 
use of Se has an additional advantage in that it can be directly expressed in terms of the 
Wigner function. 

Consider an arbitrary time dependence of the system frequency u(t). This will impose a 
constant change on the basis of coherent states (which are defined with respect to the Fock 
space of an oscillator with instantaneous frequency 00 (t)). It is therefore convenient to use 
the preferred basis (e.g., the set of coherent states in the Fock space with w(0)). The master 
equation for the Wigner function W(a, a*; t), with a, a* always referring to this basis, reads 
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32] 



^W(a, a*, t) = 2u(t)m U*^) W(a, a*, t) + ^^Re U^-\ W(a, a*, t), (15) 



or 



^W(x, y, t) = u(t) - y^-j W(x, y, t) + ^ ln ^ (t) ( x __ _ ) W (x, y, t). (16) 



/ d__ d_\ 

\ dx ^ dy J 



We omitted the diffusion terms, which describe decoherence (including relaxation). 
Therefore these equations are valid for the adiabatic stages of the cycle (assuming that 
the system has no inrinsic sources of decoherence). Equation (f!6|) is a first-order lin- 
ear equation and can be solved by the method of characteristics: using the Ansatz 
W(x,y,t) = W(x(t),y(t)) we find from ( Tl6l) the characteristic equations, 



ux — z—y- 



(17) 



dx uj dy 
— = — x — uy; — 
dt 2lu dt 2tu 

The evolution of the Wigner function W(t) = W(x (x,y,t),y (x,y,t)) is due to the initial 
distribution being "dragged" along the characteristic curves. The total phase volume occu- 
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32|, 



pied by the system obviously will not change during such an evolution. As shown in 
in the limit of fast frequency change, when the tu-terms dominate, these equations lead to 
the squeezing of the Wigner function, while in the quasistatic limit they simply describe its 
rotation as a whole, without disrupting an equilibrium state. For example, an instantaneous 
change of the oscillator frequency u c — > u>h will transform the thermal state into a squeezed 



thermal state (see 



371 ]. Eq.(4.13)) characterized by the squeezing parameter s = 0Jhjuj c 



W 



ST, 



x, y ) = —3 exp 
7rn 



— — (x 2 s + y 2 1 s 
n v 



(18) 



lc(h) 



W(x,y,t), (19) 



On the isothermic stages, the r.h.s. of Eq. ( TIB"]) must also include the diffusive term 31 1 

4 \dx 2 dy 2 J dx dy 
where r y c (h) 1S the decay rate of the corresponding reservoir. Then the equations cannot be 
solved analytically and must be dealt with numerically. 

The energy of the system and the squeezing coefficient are expressed through the Wigner 
function as 

E = (H) =coJJ dxdy (x 2 + y 2 ) W{x, y, t); (20) 

{ ff dxdy (xcos6 + ysm6) 2 W(x,y,t)\ 
s(t) = max, {„(*)} = max, | ;/ ^ — _ ^— ^ j . (21) 

In the last expression we take into account that the "cigar" of the squeezed state rotates in 
the phase plane xy with frequency u(t). 

In order to calculate the energy entropy (}Tlj) we use the expression for the Wigner function 
of a Fock state (e.g., Eq.(4.4.91) |m|): 

\n)(n\ <-> W^(a,a*) = ^^ e - 2|a|2 L n (4|«| 2 ). (22) 

7T 

Here L n (x) is the Laguerre polynomial. Finding the diagonal part of the density matrix 
in energy representation is now straightforward. It can only depend on the angle-averaged 
Wigner function 

W(\a\ 2 ) = \^ dQ W(\a\e ie , \a\e~ ie ). (23) 

27T JO 

Expanding a well-behaved function of variable r\ = 4|a| 2 , f(rf) = W(tj/4) exp[r//2], in the 
Laguerre series, 

°° POO 

f( V ) = £ A n L n { V ), A n = d V f( v ) L n ( V ) e~\ (24) 

n=0 J ° 

we see, that 

W(\a\ 2 ) = ±^p^A n W n F (\a\ 2 ). (25) 

Therefore, the diagonal matrix elements of the density matrix in the energy representation 
are p n = ir(—l) n A n /2. Writing 

dada* = \a\d\a\d6 = -dr)d9 (26) 

8 

and using (f24l) and (f23l) . we obtain 

p n = 2(-l) n J dada* e- 2|a|2 L n (4|a| 2 ) W{a, a*). (27) 

8 



In terms of x = Re a, y = Imoe, this is 



p n = 2(-l)" / / dxdye- 2 ^ 2+ ^ L„(4(x 2 + y 2 )) W(x,y). (28) 

These expressions allow a direct calculation of the diagonal elements of the density matrix, 
{p n .}n=o..oo, and with that, of S E [p\- 

One can also introduce the quasiclassical entropy, S qc [p;x\, utilizing the fact that the 
Wigner function W(p, q) reduces to a classical probability distribution after being averaged 
over the scale of ApAq > h (in our case, AxAy > 1) c 



S qc [p; x} = - J J dxdy y J J dx'dy'W(x', y')x(x -x',y- y') 

dx'dy'W(x', y')x(x -x',y- y') 



x 



In 



(29) 



Here x{ x i u) is a normalized sampling function peaked at zero, with a support dimension 
exceeding unity (e.g., x — (vrd 2 ) -1 exp[— (x 2 + y 2 )/d 2 }, d > 1/2). 

As a relevant example, let us calculate the energy entropy of the squeezed thermal state. 
From Eq.( fl8]) we can write 

,2k r 27T ,„ r°° 



Pn = 2(-l) n — [ n d6 [ drre- 2r2 - 2r2 ^ scos2e+ ^^ s[n2 ^L n (4r 2 ), 

71 JO JO 



(30) 



which can be rewritten as 



Pr. 



-l) n 



2tt 



2n roo 

d6 I dr/exp 
Jo 



H \-k(s-1/s) cos 2 9 

s 



}L n (v). (31) 



Using the formula 



drjL n (r])e- Qr > 



(Q-i) 
Q n+l 



-l)^Q(-l)^- 1 , 

9=0 



we see, that 



« f 27T ,n {|[l-7-^-l/^)cOS 2 ^]} 
Pn = / d6— 

2nJo [\ [l + f + k(s - 1/s) COS 2 } 



n+1 ' 



or 



k r 2lT n 
2tt Jo q=0 



H \-k(s-1/s) cos 2 9 

s 



q—n—l 



(32) 



(33) 



(34) 



Using the table integral (for c, d > 0) 



p2tt 

Jo 



d9 



2tt 



(c + dcos 2 9) m c m 



2 



2 Fi ( -,m; 1; , , 



we finally obtain from Eq. ([3 



Vr, 



(35) 



and the corresponding expression for Se- 

To check this, consider the non-squeezed case, s = 1. Then the hypergeometric function 
in ( 1351) is 2 Fi (...,...;...; 0) = 1, and the expression reduces to 

Ir ^ |2 - (1+K ,r^(i^)\ 

Substituting k = tanh(o;/2T), we see that indeed the populations reduce to their equilibrium 
values, 

p n = (1 - e-"' T )e-™l T = (36) 

The expression ( l35i) . while exact, is not very illuminating. A useful approximation for 
small squeezing [k (s — 1) -C 1] can be obtained directly from Eq. ([33l) rewritten as 



Vr, 



2k 



K+l [1 + K 



n+1 • 



2nJo {l- T ^;[l-l/s-(s-l/s)cos 2 ^ 

The expression in brackets is simply p^ q , the equilibrium population given by Eq. ( 136]) . 
For small values of k(s — 1), powers in the integral can be replaced by exponents, 



(37) 



1 r 2v 

Pn~p C n q — / d0exp 
27T JO 



n- 



k(s- 1) 

1 - K 



1 - 2 cos e 



exp 



n + 1 A s !) fi_2co8 2 g 
y 1 —1— /<■ V 



1 + K 



(38) 



Since 

we eventually find 

p n w p* q / k(s - V, 



dx e z 



ttIq(z) = irl {-z), 



n n + 1 
1 + 1 

1 — K 1 + K 



uj(1 - s) 
T 



(n+1/2) 



(39) 



(the last simplification works for w/T <C 1). Either approximation in ( 139]) satisfies the 
normalization condition: J2 n Pn ~ 1- 

Squeezing depopulates low-energy states and populates the high-energy ones. For small 
squeezing, the change in energy is 



8E = y~] n co 5p n ~ (s — 1) 



(40) 
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FIG. 2: (Color online.) Energy entropy Se = — YnPn^Pn °f a squeezed thermal state as a 
function of the squeezing parameter s for the ratio uj/T = 0.5. The populations of Fock states are 
given by the exact expressions ([35]) (red, lower curve) and the approximation (f39j) (black, upper 
curve) . 

[since Iq{z) = 1 + z 2 /4 + ...]. The effect on the energy entropy becomes 

SS E ~ " E( S Pn) InPn ~ " £(<*p n ) ln Pn ~ (* ~ if- (41) 
n n 

This is actually the case, as can be seen in Fig. EJ This increase in the energy and the 
entropy production (and the resulting increase in the heat transfer to the reservoirs) will 
reduce the efficiency of the Otto cycle compared to the quasistatic case. 




III. SIMULATIONS 

The equation for the Wigner function W(t) has both a drift term (|T6l) and a diffusion 
term ffl9|) and formally coincides with the Fokker-Planck equation 

dP/dt = V [PVU(x, y, t)] - D(t) V 2 P (42) 

where P is the probability density of a Brownian particle to be near the point (x, y), U(x, y, t) 
is the time-dependent potential describing deterministic forces acting on the particle, while 
D is the diffusion coefficient. It is well known that the Fokker-Planck equation is equivalent 
to the stochastic Langeyin equation (x, y) = —VU + £, with stochastic force £ satisfying the 



following conditions 



39|: 
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FIG. 3: (Color online.) Energy E of the system as a function of the reduced time, t/T^, during 
the heat engine cycle (\R\ \ < I^D- 

Therefore, the equation for the Wigner function W(t) should be equivalent to a set of 
stochastic equations: 

u\ u 7c(/i)(*) , t 
x = oj(t)y x — — x + f x 



y = -u{t)x + —y ^— + £ v , 

where the effective diffusion constant is defined as 

D = (l/4) 7c(fc) (t) (l + 2/ {exp [u{t)/T c(h) {t)\ - l} 



(43) 



(44) 



Using the Langevin equations ( 1431) we can easily simulate the evolution of the Wigner 
function. Starting at t — from the thermal Wigner function (I18p. we randomly spread 

Brownian particles according to this Gaussian probability density. Then we numerically 
monitor the evolution of each of these Brownian particles with increasing t. The obtained 
distributions of particles at any t allow us to numerically estimate the time dependence of 
the Wigner function and calculate the energy of the system E = u(t)(x 2 + y 2 ) and the 
squeezing coefficient /3 = (x 2 )/ (y 2 ), when arbitrarily changing u with time. The results are 
shown in Figs. El HI 

The energy entropy of the system, S^(t), is calculated substituting (128]) in Eq. ffT4"j) . We 
also calculate the quasiclassical entropy, S qc (t), from Eq. (l29|) . In the latter case instead 
of using the sampling function x we directly counted the number of representing particles 
within cells [xi + Ax,yi + Ay] (i.e., using a grid). The results are shown in Fig. |5l 
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FIG. 4: (Color online.) (Top) Energy E of the system as a function of time during the heat 
pump cycle \R\\ > \R2\- As the rate of adiabatic expansion decreases, the system approaches 
the operation of a quasistatic Otto cycle (inset). (Middle) The ratio ({x 2 )/(y 2 }) of quadrature 
dispersion rates during the cycle, for the same expansion rates as above. Oscillations are due to 
the rotation of the Wigner function in the phase plane with the instantaneous eigenfrequency of 
the oscillator. (Bottom) Wigner function in the phase plane at different points (A, B, C, D) of the 
Otto cycle (also shown in Fig. 1). 



The squeezing due to the finite expansion/compression rate is the only source of ineffi- 
ciency in our model system. Compared to the quasistatic case, it leads to increased entropy 
and energy at the points B and D. Thus, e.g., in the heat engine regime, decreasing the net 
work performed by the system and increasing the net heat transfer to the cold reservoir. 

From Fig. H] we see that the performance of the quantum Otto engine fast approaches 
its quasistatic limit as the squeezing of the quantum state of the working body (tunable 
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FIG. 5: (Color online.) Energy (a) and entropy (b,c) of the quantum Otto engine during the cycle 
shown in Fig. 1(b). Inset: time dependence of the oscillator frequency, to. 



oscillator) is decreased. This is consistent with the effects of squeezing being quadratic in 
the squeezing parameter, though the roughness of calculations does not allow to confirm the 
exact functional shape of this dependence. 



IV. CONCLUSIONS 



We have analyzed a simple model of a quantum Otto engine, which can be realized based 
on, e.g., Josephson devices. The roles of the working body and hot and cold reservoirs are 
played by oscillators (with tunable or fixed frequencies, respectively). We have shown that, 
depending on the relation between temperature and resonant frequencies of the reservoirs, 
the system can work either as a heat engine or a heat pump. Using the method of Wigner 
functions, we found that the source of inefficiency of this device is in the squeezing of the 
quantum state of the working-body oscillator. In particular, we found an explicit expression 
for the energy entropy of a squeezed thermal state. Though inevitable for any finite speed 
of operation, the effect of small squeezing s ~ 1 is only proportional to (1 — s) 2 . 
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Comparing the efficiencies of the quantum and thermal engines, one can say that a 
quantum Otto engine may operate with a high efficiency even at low temperatures, unlike 
the classical case. The losses can be minimized by reducing the squeezing of the working 
body. It may be worth considering incorporating quantum Otto engines in superconducting 
qubit registers as additional coolers. 
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